Binding of Nucleobases with Single- Walled Carbon Nanotubes 
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We have calculated the binding energy of various nucleobases (guanine (G), adenine (A), thymine 
(T) and cytosine (C)) with (5,5) single-walled carbon nanotubes (SWNTs) using ab-initio Hartre- 
Fock method (HF) together with force field calculations. The gas phase binding energies follow the 
sequence G > A > T > C. We show that main contribution to binding energy comes from van-der 
Wall (vdW) interaction between nanotube and nucleobases. We compare these results with the 
interaction of nucleobases with graphene. We show that the binding energy of bases with SWNTs 
is much lower than the graphene but the sequence remains same. When we include the effect of 
solvation energy (Poisson-Boltzman (PB) solver at HF level), the binding energy follow the sequence 
G>T>A>C>, which explains the experiment! 16 that oligonucleotides made of thymine bases 
are more effective in dispersing the SWNT in aqueous solution as compared to poly (A) and poly 
(C). We also demonstrate experimentally that there is differential binding affinity of nucleobases 
with the single-walled carbon nanotubes (SWNTs) by directly measuring the binding strength using 
isothermal titration (micro) calorimetry. The binding sequence of the nucleobases varies as thymine 
(T) > adenine (A) > cytosine (C), in agreement with our calculation. 



PACS numbers: 

INTRODUCTION 

Single-walled carbon nanotubes (SWNTs) are one- 
dimensional systems with different diameters and chi- 
ralities. They have large surface areas [TJ [H |3J with 
the electrons in intimate contact with the environment. 
The superb mechanical |3] and electrical [5] properties of 
SWNT have potential for many applications such as na- 
noelectronics [H [7] , actuators j8] , chemical [9j [10] , and 
flow [TT] sensors. Recently, there are potential appli- 
cations of DNA coated CNTs as a biosensor [T5J [T5] . 
Single-stranded DNA (ss-DNA) coated SWNTs field ef- 
fect transistor (FET) has been used to detect various 
odors [10], DNA hybridization [14] and conformation 
changes in DNA |15j in presence of ionic concentration. 
It has been also reported that DNA can be inserted into 
CNTs, which can be further potential applications of this 
nano-bio system. 

During synthesis SWNTs appear in bundles, but for 
their potential application we need isolated SWNT. It 
is, therefore, a big challenge to disperse SWNT bundles 
in aqueous solution. It has been reported that single- 
walled carbon nanotube bundles are effectively dispersed 
in water on sonication in the presence of single stranded 
DNA (ssDNA) [T5]. The dispersity is dependent on the 
oligonucleotide sequence. Jagota et al. [TB] has shown 
that poly (T) has higher efficiency in dispersion the nan- 
otubes compared to ploy (A) and ploy(C). Therefore, the 
wraping of carbon nanotubes by ss-DNA is sequence de- 
pendent. It was also found that the best separation 
was obtained with sequence of repeats alternating G and 



T (poly (GT)) [T7]- Studies on the nature of interac- 
tion of ssDNA, dsDNA and oligonucleotides with nan- 
otubes have made use of classical molecular dynamics 
simulations, besides experimental techniques such as ion 
exchange chromatography (IEC) |17j . atomic force mi- 
croscopy (AFM), resonance Raman spectroscopy (RRS), 
photoluminescence (PL) [TH], linear dichroism (LD) [T3] 
and directional optical absorbance [20]. There has, how- 
ever, been no estimation of the strength of interaction 
of the individual nucleobases, adenine (A), cytosine (C), 
thymine (T) and guanine (G) with the SWNTs. The ob- 
jective of this study is to determine the binding strength 
of the different nucleobases with SWNTs using ab-initio 
quantum chemical as well as classical force field calcu- 
lations. We have also determined the binding strength 
of thymine and cytosine with the SWNTs by employ- 
ing isothermal titration calorimetry (ITC), which has 
emerged as a powerful tool to study the thermodynam- 
ics of protein-protein [21] , DNA-protein [22] and protein- 
lipid interactions [23J . The advantage of ITC over kinetic 
methods is that the same experiment provides the bind- 
ing constant as well as the enthalpy of binding. In our 
study we have titrated aqueous SWNT solutions against 
aqueous solutions of nucleobases. Our results show that 
the binding affinity of thymine with SWNTs is higher 
than adenine and cytosine. We believe that the present 
study, besides providing valuable insight into the inter- 
actions of SWNTs with nucleobases, may be of potential 
use in nanotechnology. 
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FIG. 1: I, II, V & VI show the cross-sectional view of op- 
timized structure of C, G, A and T nucleobases bonded to 
nanotube. Ill, IV, VII & VIII show the lateral view of the 
corresponding optimized structures. 



THEORETICAL CALCULATIONS 

We consider a (5, 5) SWNT containing 5 unit cells 
in presence of different nucleobases (C, G, A, T) as in- 
put configuration for our quantum chemical calculations. 
The bases are kept parallel to the nanotube axis. Dur- 
ing the optimization process, all the atoms are free to 
relax. The optimized structures of nucleobases, SWNT 
and combined system are obtained using Jaguar compu- 
tational package [21] • The binding energy of nucleobases 
with SWNT is evaluated using ab-initio restricted open 
shell Hartree-Fock (ROHF). In all the calculations, we 
have used double-zeta basis set with polarization func- 
tion (6-31g**). The optimized geometries for the com- 
bined systems are shown in Fig. 1. We notice that a 
particular orientation of the base with the nanotube, as 
shown in Fig. I is achieved irrespective of different input 
configurations. Hartree-Fock (HF) or density- functional 
theory (DFT) describes the charge transfer and chemical 
bonds between the molecules but they are inadequate to 
describe the van der Waals (vdW) interaction [25J, [26] . It 
is known that the vdW is the main interaction between 
organic molecules and inert surfaces like graphite [23112"?] . 
Here we have calculated the vdW energy between nan- 
otube and different bases using classical MSCFF force 
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FIG. 2: The open circle with line shows the DOS for the com- 
bined system (SWNTmucleobase) , whereas the solid line and 
stars with dashed line show the DOS for pristine nanotube 
and nucleobases, respectively. 



field [2B] as well as AMBER 7 29J. The solvation en- 
ergies [3H1 EI] were calculated using Poisson-Boltzman 
(PB) solver of Jaguar which takes into account the pres- 
ence of continuum dielectric medium like water. The 
total binding energy is thus obtained by adding the HF, 
vdW and solvation energy. We show that the binding 
energy between the nanotube and DNA bases is mainly 
governed by vdW interaction, whereas the sequence of 
binding is influenced by the solvation energy. 

Ab-initio Calculations 

To investigate the binding energy of different bases 
with nanotube, we have calculated the energy of op- 
timized structures of the nanotube, the bases and the 
combined system separately using Hartree-Fock (HF-6- 
31g**). The binding energy of a base is calculated by 
subtracting the energy of an isolated nanotube and iso- 
lated base from the energy of the combined system i.e 

E binding — B adsorbed system ~ (E isolated nanotube E 

isolated base)- In this scenario the binding sequence is C 
> G > A > T and the corresponding binding energies 
are -0.07, -0.04, -0.03 and -0.02 eV, respectively. 

To investigate the origin of binding of nucleobases with 
nanotube we will examine the different energetic contri- 
butions, Mulliken charge analysis and electronic density 
of states (DOS). It is seen that there is no noticeable 
contribution from covalent or Coulombic energy to the 
attraction between nanotube and bases. The only no- 
ticeable contribution in the binding energy obtained from 
HF (6-31g**) comes from the exchange energy term. Nu- 
clear repulsion (E a), total one electron term (E e) and 
total two electron term (E i) are the three components in 
binding energy obtained from HF (6-3Ig**), where the 
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total two electron term consists of Coulomb (E j) and 
exchange energy (E x)- The corresponding energies for 
cytosine bonded nanotube are E a = 62595.74 eV, Ee = 
-125302.13 eV and E , = 62706.32 eV with total energy 
E binding = E A + E E + E| = -0.07 eV. As HF does not 
give the exchange energy separately we have calculated 
the exchange energy with correlation term (E xc) for cy- 
tosine bonded nanotube using density functional theory 
(DFT) at the theoretical level of B3LYP/6-31g**. The 
corresponding energies are E a = 68743.93 eV, E e = 
-137595.49 eV, E j = 68851.62 eV and E X c = -0.13 eV. 
Therefore E| = Ej+E X c = 68851.49 eV and the total 
energy E binding = Ea + Ee + E|= -0.07 eV, which is 
same as the binding energy obtained from the HF. The 
binding of cytosine with nanotube can be explained only 
if we consider the E xc term (-0.13 eV), because the sum 
of 1st three components of energies (E a + E e + E j) is 
a positive contribution of 0.06 eV. The sum E a + E e 
+ E i energy would be positive if there is no strong cova- 
lent bond between nanotube and cytosine. The positive 
value arises from purely repulsive Pauli barrier. It will be 
explained in the next paragraph that indeed there is no 
noticeable covalent bond between the nanotube and nu- 
cleobases. Therefore, the binding energy calculated from 
HF suggests that the weak attraction between nanotube 
and nucleobases comes mainly from the exchange energy. 

The DOS of pristine nanotube, isolated bases and the 
combined system are shown in Fig. 2 which clearly shows 
that there is no shift in the DOS of combined system with 
respect to the isolated SWNT and the DOS of combined 
system is the sum of DOS of pristine nanotube and the 
corresponding base. We note that the (5, 5) nanotube 
has an energy gap of 4 eV due to the finite length of 1.1 
nm. The nearest distances in optimized structures (see 
Fig. 1) between the nanotube to different nucleobases are 
C-0 3.2 A (cytosine and guanine), C-H 3.3 A (adenine) 
and C-H 3.4 A (thymine), which certainly exceed the 
sum of the covalent radii of C-0 1.5 A and C-H 1.1 A. 
Therefore, the bond between the nanotube and the bases 
is unlikely to be covalent. 

The adsorption energies in Hartree-Fock calculations 
are of similar magnitude as calculated for adenine ad- 
sorbed on Cu (110) [32], where the mutual polarization 
and Coulombic interaction between the molecule and Cu 
substrate determines the binding energy. Guided by this, 
we have carried out the Mulliken charge population anal- 
ysis to evaluate the charge transfer and charge distribu- 
tion in the combined system. As shown in Fig. 3, the ad- 
sorption does cause a minor charge redistribution. How- 
ever, there is no charge transfer between the nanotube 
and nucleobases. Out of all the four bases, the most pro- 
nounced case is that of cytosine bonded nanotube where 
the maximum charge redistribution of the order of Aq = 
+ O.Ole occurs on a carbon atom of nanotube in close 
proximity to the nitrogen atom of the cytosine (see Fig. 
3). This results an overall dipole moment in nanotube 
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FIG. 3: The open star with line shows the charge distribution 
on individual carbon atoms along the nanotube axis near (a) 
the nitrogen of C, (b) hydrogen of G, (c) nitrogen of A, (d) 
oxygen of T. The cross with line shows the averaged charge on 
each carbon atom in a ring along the nanotube axis for neu- 
tral nanotube. The open circle with line shows the averaged 
charge on each carbon atom in a ring along the nanotube axis 
for nanotube bonded with different nucleobases. 



of |P| = 0.03 Debye (D = 3.33 x 10~ 30 Coulomb me- 
ters). Since this charge redistribution is much smaller in 
magnitude compared to the adenine on copper [32], we 
can exclude any noticeable contribution to binding en- 
ergy due to mutual polarization between the nanotube 
and the bases. However, this mutual polarization can 
explain the origin for the specific orientation of the bases 
with respect to the nanotube in HF optimized structures. 
The optimized structures (Fig. 1) indicate that cyto- 
sine binds to nanotube surface via nitrogen atom of the 
pyrimidine ring with partial contribution from the oxy- 
gen of C=0 group. Similarly for guanine, adenine and 
thymine, the binding occurs via either the nitrogen atom 
of the pyrimidine ring or the oxygen atom of the carbonyl 
group. As these nitrogen and oxygen atoms have more 
electron affinity, the resultant dipole moment of the pure 
bases is determined by these atoms. 

Classical Force Field Calculations 

Since the vdW interaction is the dominant interac- 
tion between organic molecules and nanotubes, we have 
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FIG. 4: (color online) Classically optimized geometry of nu- 
cleobases on nanotube: (a) cytosine, (b) guanine, (c) adenine 
and (d) thymine. 



carried out classical force field calculation to calculate 
the vdW energy between the nanotube and nucleobases. 
Here we have used MSCFF force field parameter [28] be- 
tween the different atoms of nanotube and nucleobases. 
The input configuration is the same as that in the quan- 
tum chemical calculations. From the optimized struc- 
tures of combined systems, we note that all the nucle- 
obases remain parallel to the nanotube at a distance of 
3.25 A as shown in Fig. 4. The parallel configuration is 
obtained because of maximum vdW interaction between 
the substrate and the molecule ans the stacking arrange- 
ment (see Fig. 4) is similar to AB stacking of graphite 
layers but it deviates frob perfect AB stacking, which 
is more prominent in G, due to five and six-membered 
rings of bases and nanotube curvature effects. The vdW 
binding energies are calculated by subtracting the vdW 
energy of isolated systems from the energy of the ad- 
sorbed system and the corresponding energies are -0.54, 
-0.51, -0.47, -0.39 eV for G, A, T and C, respectively 
and the vdW energy calculated using AMBER 7 remains 
similar to that of MSCFF (see table 1) . Now the resul- 
tant gas phase energies are obtained by adding the HF 
and vdW energies and the sequence becomes G > A > 
T > C with energies -0.58, -0.54, -0.49, -0.46 eV, respec- 
tively (see table 1). It is clear that vdW energy has much 
more contribution to binding than the HF exchange en- 
ergy. We note that the gas phase binding sequence of 
nucleobases with SWNT remains same as the binding se- 
quence of nucleobases with graphite [33] . 

To compare the adsorption of nucleobases to nanotube 
and graphite, we have calculated the binding energy of 
all the bases with graphite. For adenine, the calculated 
binding energies using HF (6-31g**) and classical force 



field (vdW) are -1.50 kcal/mol (-0.06 eV) and -18.69 
kcal/mol (-0.8 eV), respectively. Therefore, the total 
binding energy (HF+vdW) of adenine on graphite is - 
20.19 kcal/mol (-0.87 eV). These energies are very similar 
to the value of -0.07 eV calculated using DFT with Gen- 
eralized Gradient Approximation (GGA) and the total 
energy of -1.07 eV (GGA + vdW) [26]. This compar- 
ison is very important as we calculate the the binding 
energy using HF and classical force field separately but 
it givesas accurate as DFT (GGA + vdW) calculation. 
The HF optimized structure for graphene-adenine com- 
bined system is shown in Fig. 5(a) and 5(b). We see that 
unlike nanotube adenine remains almost in parallel with 
graphene. The vdW interaction for graphite with other 
bases G, T and C bases are, respectively, -0.83, -0.75 and 
-0.7, eV. These values compare well with the reported 
values [33J [M] of graphite-base interaction using molec- 
ular dynamics and DFT calculation. We note that vdW 
interaction for the nanotube is smaller than the graphite 
due to curvature effects. The clasically optimized struc- 
tures for different bases with graphene are shown in Fig. 
5(c), (d), (e) and (f) for C, G, A and T respectively. Like 
in nanotube, 7r-orbitals of nucleobases and graphene min- 
imize their overlap to lower the repulsion interaction giv- 
ing rise to AB stacking like arrangement, which is more 
perfect than nanotube. 



Solvation Energy Calculations 

The binding energies calculated so far are in gas phase 
condition. In real experimental situations, the nucle- 
obases bind with the nanotubes in aqueous solution like 
water. To calculate the binding energy in solution phase 
we have carried out solvation energy calculations using 
Poisson-Boltzman (PB) [30] solver of Jaguar package at 
the HF/6-31g** level, where the gas phase optimized 
structures (HF) are taken as input and the solvation en- 
ergies (S.E) are calculated without optimizing the geom- 
etry further in presence of water. Here we have assigned 
the van der Waal radii for different atoms according to 
the reference [35]. Now the contribution of solvation en- 
ergy to the binding energy is calculated by subtracting 
the solvation energy of isolated systems from the solva- 
tion energy of adsorbed system and the energies are -0.02, 
0.01, 0.05 and 0.16 eV for T, G, A and C, respectively 
(see table 2). The total binding energy is obtained by 
adding the HF, vdW and the solvation energy. The rela- 
tive binding for the monomeric bases isG>T>A>C 
with energies -0.57, -0.51, -0.49 and -0.3 eV, respectively. 



EXPERIMENTAL DETAILS 

SWNT bundles were prepared by the arc discharge 
method followed by purification process, involving acid 
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FIG. 5: (color online) (a) HF optimized structure for 
graphene-adenine combined system, (b) cros-sectional view. 
Classically optimized structures for (c) cytosine, (d) guanine, 
(e) adenine and (f) thymine. 



washing and high-temperature hydrogen treatment |36j . 
The nanotubes are characterized by the thermogravimet- 
ric analysis, transmission electron microscopy (TEM), 
Raman spectroscopy [RS] and near infrared spectroscopy 
(NIR) [36\ . The average diameter of the nanotubes was 
about 1.4 nm, the length being a few microns. 100 /ig of 
the SWNT sample was suspended in 1.6 ml double dis- 
tilled water having a pH value of 6.7. The solution was 
sonicated for 5 hrs at a power level of 3W. After the son- 
ication, the solution was kept for a day to sediment all 
the big bundles of SWNT. The supernatant solution was 
stable for more than 15 days. The nucleobases obtained 



TABLE I: Theoretically calculated values of interaction ener- 
gies of monomeric nucleobases with SWNT 



Nucleo 
base 


HF 

eV 


vdW(eV) 
MSCFF 


HF+vdW 
eV 


S.E 
eV 


Total 
eV 


vdW(eV) 
AMBER 


Guanine 


-0.04 


-0.54 


-0.58 


0.01 


-0.57 


-0.6 


Thymine 


-0.02 


-0.47 


-0.49 


-0.02 


-0.51 


-0.5 


Adenine 


-0.03 


-0.51 


-0.54 


0.05 


-0.49 


-0.53 


Cytosine 


-0.07 


-0.39 


-0.46 


0.16 


-0.3 


-0.45 



from the Sigma Aldrich Chemicals were used as received. 
The binding of nucleobases with SWNTs was achieved by 
mixing of 300 fil of 10 mM nucleobases solution (prepared 
in double distilled water with pH of 6.7) with 1.6 ml of the 
above SWNT solution. We have chosen to suspend the 
nanotubes and the nucleobase in water instead of a buffer 
solution because the latter may influence the surface of 
the nanotubes. We could not perform the experiment 
with guanine in aqueous solution due to its insolubility 
in water. Calorimetric titrations were performed at 5 °C 
with a Micro-Cal VP-ITC instrument. We checked the 
stabilization of the instrument by performing a water- 
water run, the heat change so obtained being less than 
0.6 fiCal per 10 /A injection of water. Here a rotating 
stirrer-syringe (270 rpm) containing 300 /il of lOmM nu- 
cleobase solution injects in equal steps of 10 /il solution 
at 180-sec intervals into a cell containing 1.4 ml of SWNT 
solution until saturation is reached. 



EXPERIMENTAL RESULTS 

The ITC response recorded during the titration of nu- 
cleobases with the SWNTs are shown in Fig. 6. Fig. 
6a, 6d and 6g show the ITC data for the blank titration 
of 10 mM thymine, cytosine and adenine, respectively 
against double distilled water. Fig. 6b, 6e and 6h show 
the raw data of heat of reaction due to nucleobase-SWNT 
binding during each injection for thymine, cytosine and 
adenine, respectively. We can sec that the blank titration 
of adenine (see Fig. 6g) with double distilled water is an 
endothcrmic reaction with higher magnitude of heat of 
reaction. The raw ITC response for adenine (see Fig. 6h) 
with SWNT solution also shows an endothermic response 
with different magnitude of heat of reaction. Therefore, 
the heat of reaction of adenine with SWNTs is obtained 
by subtracting the blank titration from the raw titration, 
as shown in Fig. 6i. The nature of the binding curve 
(Fig. 6b, 6e and 6i) arises due to lowering of accessible 
surface area of SWNT in successive injections as nucle- 
obases bind with the SWNTs. The exothermic tail in the 
main curve is due to the dilution effect of nucleobases 
during the titration as there is no free surface of SWNTs 
available to bind. We have also seen that the heat of 
dilution of SWNT solution with water is less than 5 cal 
per 10 /il injection of water. The raw ITC of thymine 
(Fig. 6b) shows the appearance of second peak in each 
injection, which may be as a result of more available free 
surfaces due to possible removal of entanglement between 
the nanotube bundles. Fig. 6c, 6f and 6j correspond to 
the integrated heat change (enthalpy change in cal/mol) 
due to each injection of thymine, cytosine and adenine, 
respectively after correcting for heat of dilution (at each 
injection) and are plotted as a function of injected nucle- 
obase volume. The experiments were repeated twice and 
show good reproducibility as seen by the relatively small 
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FIG. 6: The blank titration for T (a), C (d) and A (g). The raw ITC response for T (b), C (e) and A (h). (i) Exothermic 
reaction of adenine with SWNTs after subtracting the blank titration from the raw titration, (c), (f) and (j) show the integrated 
heat of reaction at each injection for T, C and A, respectively. 



standard errors. 

The strength of interaction of the nuclcobases with 
SWNTs is evaluted from the exothermicity (integrated 
heat in Fig. 6c, 6f and 6j) of the first injection of nucle- 
obases with SWNTs and the binding sequence is T (-320 
cal/mol) > A (-219 cal/mol) > C (-134 cal/mol). Since 
the heat of reaction curve in Fig. 6 are not in sigmoidal 
nature, the absolute value of binding energy can not be 
calculated from Fig. 6. This is because of the limited 
solubility of the nanotubes as well as its bundle effect in 
aqueous solution. 



CONCLUSIONS 

Our results from HF combined with forcefield calcu- 
lation indeed suggests that the binding of nucleobases 
with SWNTs is mainly due to weak vdW force and the 
gas phase binding sequence isG>A>T>C like 
the binding sequence of nuclcobases with graphene. We 
also show that compared to graphene, nanotube has 
lesser amount of binding affinity with nucleobases duo 
to curvature effects. We have shown that the solution 
phase binding sequence of nucleobases with nanotube is 
G > T > A > C, which explains the earlier observations 
that oligonucleotides made of thymine bases are more 
effective in dispersing the SWNT in aqueous solution 
as compared to poly (A) and poly (C). We have shown 
by direct calorimetry that there is differential binding 
affinity in the interaction of nucleobases with SWNTs, 
the sequence being T > A > C, in agreement with the 



calculation. 
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